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Abstract 


The inclusive transverse momentum and rapidity distributions for top 
quark production at the Fermilab Tevatron are presented both in order 
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all orders of QCD perturbation theory. The resummed results are uniformly 
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I. INTRODUCTION 


At present the top quark has not yet been conclusively discovered at the Fermilab Teva- 
tron, even though there are events which look similar to those expected from top quark 
decays. The experimental situation is presently confusing because both the CDF [Q and DO 
1^ collaborations only have limited statistics. As more events are collected one expects that 
the situation will be clarihed. 

At the Tevatron, the top quark should be mainly produced through ti pair production 
from the light mass quarks and gluons in the colliding proton and antiproton. Both the 
top quark and the top antiquark then decay to (hF, b) pairs, and each W boson can decay 
either hadronically or leptonically. The 6-quark becomes an on-mass-shell S-hadron which 
subsequently decays into leptons and (charmed) hadrons. A large effort is being made to re¬ 
construct the top quark mass from the measured particles in the decay, which is complicated 
by the fact that the neutrinos are never detected. Also there are additional jets so it is not 
clear which ones to choose to recombine The best channel for this mass reconstruction 
is where both W bosons decay leptonically, one to a (e, Ue) pair, the other to a (/x, pair 
(a dilepton event) because the backgrounds in this channel are small. When only a single 
lepton is detected then it is necessary to identify the 6 quark in the decay to remove large 
backgrounds from the production of hF-|- jets |^. In all cases the reconstruction of the parti¬ 
cles in the hnal state involves both the details of the production of the top quark-antiquark 
pair as well as the knowledge of their fragmentation and decay products. The CDF col¬ 
laboration |l| have reported two events with dilepton hnal states, six events with a single 
lepton and a 6-quark identihed by a secondary vertex, and seven single lepton events with 
the 6-quark identihed by a semileptonic decay. The CDF collaboration then constructed 
a likelihood function for the invariant mass and quoted the value rritop = 174 ± 
GeV^c^. The top quark cross section quoted by the CDF collaboration is 13.pb. The 
DO collaboration p[ have reported on nine events with an expected background of 3.8 ± 0.9. 
If the excess is due to ti production and if the top quark mass is 180 GeV/c^, then the top 
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quark cross section is 8.2 ± 5.1 pb. 

The values for the top quark production cross section as a function of the top quark 
mass used by the CDF and DO collaborations contain both the NLO QCD corrections 
and an extension to include the resummation of initial state soft partons to all orders in 
perturbation theory |^. A recent summary of the theoretical predictions has been presented 
in |l^. In the Laenen et al. paper [§ the DIS factorization scheme was used with the 


MRSD_ parton distributions the two-loop running coupling constant with five active 
flavors, and Aqcd = 0.152 GeV. 

In the analysis of the decay distributions one needs knowledge of the inclusive differ¬ 
ential distributions of the heavy quarks in transverse momentum pt and rapidity y. These 
distributions are known in NLO ||l2| , p!3[| . What we would like to discuss in this paper is an 
updating of the resummation effects on the inclusive transverse momentum distribution of 
the top quark. In the original paper 


it was not known which mass to choose whereas 
now we can assume that the mass is 175 GeV/c^. Also we discuss here how the resummation 
effects modify the rapidity distribution of the top quark. Since there have been suggestions 
of using the mass and angular distributions in top quark production to look for physics 
beyond the standard model []^ it is very important to know the normal QGD predictions 
for these quantities. 

We hrst summarize what is known on the top quark cross section. If the top quark mass 
is 175 GeV/c^ then the dominant production channel is q + q ^ t + i. In lowest order QGD 
perturbation theory it contributes about 90 % of the total cross section with the reaction 
g + g ^ t + i making up the remaining 10%. One notes that the NLO corrections in the qq 
channel are small, whereas those in the gg channel are more than 80% . At this large top 
quark mass the qg and qg channels give negligible contributions so we do not consider them. 
Even though the gg channel contribution is small in Born approximation it can be signihcant 
in NLO due to multiple soft parton radiation. These large corrections are predominantly 
from the threshold region for heavy quark production. It was shown previously []^ that 
initial state gluon bremsstrahlung (ISGB) is responsible for the large corrections at NLO 
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near threshold. 


In the dominant logarithms from ISGB, which are the cause of the large corrections 
near threshold, were carefully examined. Such logarithms have been studied previously in 
Drell-Yan (DY) production at hxed target energies (again near threshold) where they 
are responsible for correspondingly large corrections. The analogy between DY and heavy 
quark production cross sections was exploited in [Q and a formula to resum the leading and 
next-to-leading logarithms in pQCD to all orders was proposed. Since the contributions due 
to these logarithms are positive (when all scales /i are set equal to the heavy quark mass m ), 
the effect of summing the higher order corrections increases the top quark production cross 
section over that predicted in O(a^). This sum, which will be identihed as cTres, depends on a 
nonperturbative parameter /xq- The reason that a new parameter has to be introduced is that 
the resummation is sensitive to the scale at which pQCD breaks down. As we approach the 
threshold region other, nonperturbative, physics plays a role (higher twist, bound states, 
etc.) indicated by a dramatic increase in ag and in the resummed cross section. This is 
commonly called the effect of the infrared renormalon or Landau pole []K|. We chose to 
simply cut off the resummation at a specihc scale po where Aqcd << ho << since it 
is not obvious how to incorporate the nonperturbative effects. Note that our resummed 
corrections diverge for small /xq but this is not physical since they should be joined smoothly 
onto some nonperturbative prescription and the total cross section will be hnite. Another 
way to make it hnite would be to avoid the infrared renormalon by a specihc continuation 
around it, i.e. the principal value resummation method |T^J^ . However, at the moment 
our total resummed corrections depend on the parameter /xq for which we can only make a 


rough estimate. See [14] for more details. 


II. SOFT GLUON APPROXIMATION TO THE INCLUSIVE DISTRIBUTIONS 

To make this paper self-contained we list some relevant formulae. The partonic processes 
under discussion will be denoted by 
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i{ki) +i{k2) Q{j)i)+Qij>2), 

where i,j = g,q,q- The kinematical variables 

s={ki + k2y , ti = {k2 - P 2 Y - , ui = {ki - P 2 Y - 


( 2 . 1 ) 


( 2 . 2 ) 


are introduced in the calculation of the corrections to the single particle inclusive differential 
distributions of the heavy (anti)quark. We do not distinguish in the text between the 
heavy quark and heavy antiquark since the distributions are essentially identical in our 
calculations. Here s is the square of the parton-parton c.m. energy and the heavy quark 
transverse momentum is given by pt = (tiUi/s — The rapidity variable is dehned by 

exp(2?/) = Ui/ti. The Born approximation differential cross sections can be expressed by 

s -—-= d[s + ti +ui)a^As,ti,ui), (2.3) 

(it\ dVj-i 

with 


= 'Kal{p^)KqqNCF 


rtf + u\ ^ 2m? 


(2.4) 


and 


agAsAi.ui) = 27ra;{p^)KggNCF Cf - Ca 


tlMll 


X 


vti Ui Am?s 

^U\ t\U\ 


h _ 

V f-i ?/i / 




(2.5) 


Here the color factors are Ca = N and Cf = {N"^ — 1)/{2N). The color average factors 
are Kgg = N~‘^ and Kgg = {N‘^ — 1)“^. The parameter p denotes the renormalization 
scale. In |]^ the inclusive cross section was examined near threshold (s ~ 4m^) where 
the contributions from the radiation of soft and collinear gluons are large. A variable S 4 = 
s + ti + ui was dehned, where ti = (/c 2 — P 2 Y — 'm? and ui = {ki — P 2 Y — m? are inelastic 
variables in the channel i{ki) + j(/c 2 ) —^ Qipi) + Q{.P 2 ) + g{k?). The variable S 4 > 0 now 
depends on the four momentum of the extra parton(s) emitted in the reaction. In the Born 


approximation there are no additional partons so S 4 = 0. In the NLO contributions were 
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examined in the soft region (where S4 —>■ 0) and it was found that the dominant contribution 
to the NLO cross section in this region had a similar form to the NLO correction in the Drell- 
Yan process. As the latter correction is known exactly in NNLO this correspondence 
was used to write the differential cross section in order as follows 

" dt, A. — = (;;;2j«(®4 - a) 


1=0 "^4 

1 , n. . A 




/ + ! 




( 2 , 6 ) 


Here a small parameter A has been introduced to allow us to distinguish between the soft 
(s 4 < A) and the hard (54 > A) regions in phase space. The quantities ai^jj?) contain 
terms involving the QCD /3-functions and color factors. The variables ti and ui were then 
mapped onto the variables S 4 and cos9, where 9 is the parton-parton c.m. scattering angle. 
After explicit integration over the angle 9, the resulting series was exponentiated by the 
introduction of the S 4 variable into the argument of the running coupling constant. 

As noted in the previous paper 0 in addition to the total cross section we can also 
derive the resummed heavy (anti)quark inclusive pt (and below the y) distributions. The 
transverse momentum pt of the heavy quark is related to our previous variables by 


ti = --{s - S4 - [(s - 54)^ - , 

ill = -i{s - S4 + [(s - S4)^ - 4 smy,'^} , 


(2.7) 


( 2 , 8 ) 


with ml = +Pt- The double differential cross section is therefore 




dti du 


dptds4 


(2,9) 


with the boundaries 


0<p^<- — , 0 <S 4 <s — 2 mt\fs. 


( 2 . 10 ) 


The O(a^) contribution to the inclusive transverse momentum distribution daij/dpl is given 
by 
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( 2 . 11 ) 


dpi 

1 ( '^4 


dsA 


X 


1 

/Ti 

-(jfAs,Si,pl) , 


^ ~ (^)^(^4 - A) + -— ln'+^ (—) 5 ( 54 )} 

'-S 4 / + 1 d 


[{s - s^y - Asm'lyA *4 

where we have inserted an extra factor of two so that / dpfda/dp^ = atot- After some algebra 
we can rewrite this result as 


d(^ii\s,p‘f) 


dp 


2k-l 

z=o 






ds4 — In^ 


S 4 


S 4 


X 


( dalf{s,S4,p‘f) _ dalf{s,0,p‘f) ^ 
I dpf f 


dpt 


+ 

with the dehnition 

ddlf{s,S4,Pt) 


l + l 


— 2mtS^A^ 

dalf{s,0,pl)i 

V 77^2 J 

dpi J 


( 2 . 12 ) 


dpt 


s[(s — 54 )^ — AsnifyA 


(^fj{s,S4,Pt) , 


(2.13) 


where da^f {s,0,p‘f)/dp'f = dafj^{s^pf)/dpi again represents the Born differential pt distri¬ 
bution. For the qq and gg subprocesses we have the explicit results 

o o. r.‘^\ 1 \ 


dpi ^ [(s - 34? - 4sm2]V2 

As - 34? -^SpA 


X 


(2.14) 


and 


daf^{s,S4,pl) ^ 2 ^ 2^r^ 1 1 


Cr-C 


ml 


[{s — S4? — 2sml 4m^ / m? 

+ ^(l- 2 

^ mi 


sml 


ml 


(2.15) 


Since the above formulae are symmetric under the interchange ti ui the heavy quark 
and heavy antiquark inclusive pt distributions are identical. Note that (2.12) is basically the 
integral of a plus distribution together with a surface term. 
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The corresponding formula to (2.12) for the rapidity y of the heavy quark is obtained by 


using 


ti = --^^-^^(1 - tanhy) , 


(2.16) 


Ui = — ^ (1 + tanh|/) . 

The double differential cross section is therefore 

y (Paij{s, Si,y) 

dti dui s — S 4 dy ds^ 


(2.17) 


(2,18) 


with the boundaries 

- ^ In (^^) < y < ^ In , 0 < S 4 < s - 2v/smcosh|/, (2.19) 

where = 1 —4m^/s. The O^a^) contribution to the inclusive rapidity distribution daij/dy 


is given by 


'll'] 2A: 1 ns—coshy 

E / ds, 

dy Jo 


X 


/ ^ — ^4 \ R / \ 


'2s^ cosh^ 1 /^ 

After some algebra we can rewrite this result as 


d(yf^\s,y) 

dy 


= a. 


X' 


r ns—coshy 1 / c. \ 

(y)i:a,(y) / 4 * 4 -In' ^ 

■ ddlf (s, S4, y) ddif (s, 0,y), 


dy 


dy 


1 + 1 


fS — cosh 7 / 

dalf{s,0,y)i 

V ^2 J 

dy J 


( 2 . 20 ) 


( 2 . 21 ) 


with the dehnition 


dalf{s,S4,y) 

dy 


S S4 . , 

2 ^ij ('®’ '^ 4 ) y) ) 


2s^ cosh^ y 


( 2 . 22 ) 



where dalj\s, 0, y)/dy = dalj\s, y)/dy again represents the Born differential y distribution. 
For the qq and gg subprocesses we have the explicit formulae 


d(Tqq{s,S4.,y) 

dy 


= 'Kal{g^)KqqNCF 


Sa 


2s^ cosh"^ y 


X 


cosm y 


54 )^ / ,2 • 1 2 \ 

cosh y + sinh yj H- 


(2.23) 


and 


d(Tf^{s,Si,y) 

dy 


AT^al{yi^)KggNCF 
Cf - Ca 


s — S4 


X 


s - 54 )^ 


2 s 2 cosh y 
X 


4s^ cosh^ y^ 

9>m?s cosh^ y / 4m^s cosh^ y 


(s - 84)2 


.(1- 


(s - 84)2 


cosh^ y + sinh^ y 

); 


(2.24) 


Since the above formulae are symmetric under the interchange ti <-»• ui the heavy quark and 
heavy antiquark inclusive ^/-distributions are identical. Also ( 2 . 21 ) is again of the form of 
a plus distribution together with a surface term. Finally, we note that the terms in (2.12) 
and ( 2 . 21 ) are all finite. 


III. RESUMMATION PROCEDURE IN PARTON-PARTON COLLISIONS 


The resummed contribution to the top quark cross section can be written as 

2 d (jjj(s, ti, Ui) 


dt\ du^ 


c(S 4 jji^ 


CT^j(s,fi,Mi), (3.1) 


where 


/ S 4 m \ 
\m^ jjL^ J 


A^as mA In^ ^ exp(-777E). (3.2) 

TT \m^ J ) 1 (1-|-7) 

The straightforward expansion of the exponential plus the change of the argument in as via 
the renormalization group equations generates the corresponding leading logarithmic terms 
written expicitly in [Q . The scheme dependent A and as in the above expression are given 
by 
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-4 = 2; as{y, /i^) = = 


47r 


in the MS scheme, and 

-4=1; as{y, y^) = as{yy^) = 


/3o ln(|/^/^/i^/A2) ’ 


dvr 


(3.3) 


(3,4) 


AMa^iVA^) ’ 

in the DIS scheme, where /?o = 11/3 Ca — 2/3n f is the lowest order coefficient of the QCD 
/3-function. The color factors Cij are dehned by Cqq = Cp and Cgg = Ca, and is the 
Euler constant. The quantity rj is given by 


h = 




(3.5) 


( 3 q \ dvr /i" 

Following the procedure in |T^ for the resummation of the order contributions to the 
Pt distribution we have 

daij{s,pi) ^ g daij\s,p‘f) 


dpi 


k =0 


dpi 


CLS^ 


/so 


X 


(iS4 
(o; 

'ij 


I dpi dpi 1 

/ s - 2mts^^^ m^\ da§\s,pi) 

1 / dpi 




(3.6) 


Note that we now have cut off the lower limit of the S 4 integration at S 4 = Sq because ag in 
(3.2) diverges as S 4 —0. This parameter sq niust satisfy the conditions 0 < Sq < s — 2mts^^‘^ 
and SqIvt? << 1. It is convenient to rewrite Sq in terms of the scale y as 

(iAC 


So 






(MS scheme); 


(3.7a) 


^ (DIS scheme). 


(3.7b) 


Here yo is a nonperturbative parameter satisfying Aq^jd << yl « y^- The derivative 
of f{s 4 ^/m?,m?/y‘^) is obtained from (3.2). It is equal to 
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(3.8) 


df{s4/m?, m?/jj?' 
dsA 


In ^ + r;} 

'■ 'TT Ul^ 'rri'^ ’ 


S 4 '■ TT 
X exp 




2^n 


>■ TT J i (1 + 77 ) 

X exp(-r/7i5), 


where we have neglected terms which are higher order in ag- 
The analogous formula for the rapidity distribution is 
daijjs^y) ^ g da\f{s,y) 


dy 


k=0 


dy 


rs- 2 ms^/^ coshy df {sA/, 771^/fi^) 

/ - 

J SQ cio4 

j dd^f{s,S4,y) _ dalf{s,0,y) ^ 

I dy dll j 


dy 


/s — coshy w? \da^^\s^y) 

^ w? ' dy 

with the conditions 0 < Sq < s — cosh 7 / and So/m? << 1. 


(3.9) 


IV. NUMERICAL RESULTS 


Following the notation in |]I^ the total hadron-hadron cross section in order cxj is 

(yP{S,m^) = ^j dr<Fy(r,p^)aff , (4.1) 

where S is the square of the hadron-hadron c.m. energy and i,j run over q, q and g. The 
parton flux <Fy(r,/i^) is dehned via 

dx 


Jt X X 


(4.2) 


and Hij is a product of the scale-dependent parton distribution functions f^{x,y^), where 


h stands for the hadron which is the source of the parton i 




(4.3) 


The mass factorization scale y is chosen to be identical with the renormalization scale in 
the running coupling constant. Since the pt distribution in hadron-hadron collisions is not 
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altered by the Lorentz transformation along the collision axis from the parton-parton c.m. 
frame, we can write an analogons formnla to (4.1) for the heavy-qnark inclnsive differential 
distribntion in p'f 


dpi 


E 


Ami / ^ 


dr <h 


da. 






(k) 


{tS, m‘^,p‘l,p^) 


dpi 


(4.4) 


with + Pt- Iii the case of the all-order resnmmed expressions the lower bonndaries 

in (4.1) and (4.4) have to be modihed according to the conditions sq < s — or 

Sq < s — (see above). Resnmming the soft glnon contribntions to all orders we 

obtain 


^res 


{S,vr?) = E / dT^i^[T,p?)ai^[TS,vr?,p^) , 

nA •>'^0 


(4.5) 


where is given in (3.24) of |T5 and 

(m -I- (m^ -I- 


To = 


S 


(4.6) 


with So = rn'^ipl /(MS scheme) or sq = m^(po/h^) (DIS scheme) [see (3.7)]. The 
all-order resnmmed differential distribntion in p^ is 


dpi 


= E/ dr<l)ij(r,/i' 

'TO 


daijirS, rn/,pi,p^) 


dpi 


(4.7) 


with daij/dpi given in (3.6) and 


To = 


{rrit + {mi + So)^^)" 


(4.8) 


The corresponding formula to (4.4) for the heavy quark inclusive differential distribution in 

y is 

daP(S,m^y) _y. A , ,. d4’(TS,m\y,y^) ^ 

Order by order in perturbation theory the heavy quark rapidity plots in the parton-parton 
c.m. frame show peaks away from |/ = 0 (see £g. 7 in []^). However, upon folding with the 
partonic densities the heavy quark rapidities in the hadron-hadron c.m. frame peak near 
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y = 0. Therefore we will assume that the plots for the resummed rapidity distribution show 
a similar feature. The all-order resummed differential distribution in Y is therefore taken to 
be 


/■ dr4.„(r.y) 


dY ' dy 

with daij/dy given in (3.9) and 

(mcosh?/ -|- {m? cosh^ y -|- 


(4.10) 


To = 


s 


(4.11) 


The hadronic heavy quark rapidity Y is related to the partonic heavy quark rapidity y by 

Y = y + -\n — . (4.12) 


We now specialize to top quark production at the Fermilab Tevatron where \fS = 1.8 
TeV and choose the top quark mass to be m = 175 GeV/c^. In the presentation of our 
results for the exact, approximate and resummed hadronic cross sections we use the MRSD_ 


parametrization for the parton distributions . Note that the hadronic results only involve 
partonic distribution functions at moderate and large a:, where there is little difference 
between the various sets of parton densities. We have used the MRSD_ set 34 as given in 
PDFLIB in the DIS scheme with the number of active light flavors nf = 5 and the QCD 
scale As = 0.1559 GeV. We have used the two-loop corrected running coupling constant 
as given by PDFLIB. Since we know the exact O(a^) result, we can make an even better 
estimate of the differential distributions by calculating the perturbation theory improved pt 
and y distributions. We dehne the improved pt distribution by 




dpt dpt 

and the improved Y distribution by 


l da 


( 1 ) 


da 


( 1 ) 


dY 


dY 


+ 


n 

1 \ exact 

dpt 

dpt 

|app 7 

(4.13) 


da]^^ 


(4.14) 

i'\ r exact 

dY 

dY 

|app 7 

and da^^ /dY 

1 exact 

are known and da^^ /dpt 

|app nnd 


da^^/dY |app are included in da'^ff/dpt and da'^ff/dY respectively. We note that here da^'^'^ 
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denotes the 0(a”+^) contribution to the differential cross section. Moreover , |exact 

denotes the exact calculated differential cross section, and |app the approximate one 

where only the leading soft gluon corrections are taken into account. 

First we present the differential pt distributions at y/S =1.8 TeV for a top quark mass 
m = 175 GeV/c^. For these plots the mass factorization scale is not everywhere equal to 
m. We chose /i = m in sq, fk{.Si/m? /p?) and but p = in the MRSD_ parton 
distribution functions and the running coupling constant as{p)- We begin with the results 
for the qq channel in the DIS scheme. In £g. 1 we show the Born term /dpi^ the first 
order exact result da^^ /dpt |exact, the first order approximation da^^/dpt |app, the second 
order approximation da'^^ /dpt |app, and the resummed result da^^/dpt for pq = 0.05 m 
and for po = 0.1 m. These are the same values for po that were used in p. As we 
decrease po the differential cross sections increase. In fig. 2 we show the exact O(a^) result 
da'^H /dpt + da^H /dpt |exact, and, for comparison, da'^^/dpt for po = 0.05 m and po = 0.1 m. 
The resummed differential cross sections were calculated with the cut S 4 > sq while no such 
cut was imposed on the phase space for the individual terms in the perturbation series. 
The improved differential distributions are uniformly above the exact O(a^) result. It is 
also evident from fig. 2 that the resummation of the soft gluon contributions to the pt 
distribution modihes the exact O^a^) result only slightly for the values of po that have been 
chosen. 

We continue with the results for the gg channel in the MS scheme. The corresponding 
plots are given in figures 3 and 4. In this case the values of po have been chosen to be 
Po = 0.2 m and po = 0.25 m to correspond to those in |P]. Note that po need not be the 
same in the qq and gg reactions because the convergence properties of the QCD perturbation 
series could be different in these channels and moreover depend on the factorization scheme. 
The first and second order corrections in the gg channel in the MS scheme are larger than 
the respective ones in the qq channel in the DIS scheme. In fact, for the range of pt values 
shown the second-order approximate correction is larger than the first-order approximation. 
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Hence, the relative difference in magnitude between the improved /dpt and the exact 
0 ( 0 ;^) results is signihcantly larger than that in the qq channel in the DIS scheme. 

We hnish our discussion of the differential pt distributions with the results of adding the 
qq and gg channels. The plots appear in hgures 5 and 6. It is evident that resummation 
produces an enhancement of the exact 0 ( 0 ;^) result, with very little change in shape. 

Now we turn to a discussion of the differential Y distributions at = 1.8 TeV for a 
top quark mass m = 175 GeV/c^. In this case we set the factorization mass scale equal 
to m everywhere. We begin with the results for the qq channel in the DIS scheme. In £g. 
7 we show the Born term da^^^/dY, the first order exact result da%^/dY [exact, the first 
order approximation da^^ /dY |app, the second order approximation da^^/dY |app, and the 
resummed result dY^^/dY for p,Q = 0.05 m and /io = 0.1 m. In fig. 8 we show the exact 
0(ag) result da^^/dY + da^^^/dY [exact, and, for comparison, da^'^/dY for /tq = 0.05 m 
and /io = 0.1 m. Again, the resummed differential cross sections were calculated with the 
cut S 4 > So while no such cut was imposed on the phase space for the individual terms 
in the perturbation series. It is also evident from fig. 8 that the resummation of the soft 
gluon contributions to the Y distribution modifies the exact 0 (q;^) result only slightly for 
the values of /io that have been chosen. 

We continue with the results for the gg channel in the MS scheme. The corresponding 
plots are given in figures 9 and 10. Here, the values of /io are /io = 0.2m and /io = 0.25m. As 
in the case of the pt distributions, the first and second order corrections in this channel are 
larger than the respective ones in the qq channel in the DIS scheme. For the range of Y values 
shown the second-order approximate correction is larger than the first-order approximation. 
Again, as in the pt distributions, the relative difference in magnitude between the improved 
dY^^/dY and the exact 0{a^) results is significantly larger than that in the qq channel in 
the DIS scheme. 

Finally, we conclude our discussion of the differential Y distributions by showing the 
results of adding the qq and gg channels. The plots appear in figures 11 and 12. Again, it is 
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evident that resummation produces a non-negligible modification of the exact 0{al) result. 
However, the shape of the distribution is unchanged. 

We have shown that the resummation of soft gluon radiation produces a small difference 
between the perturbation improved distributions in pt and Y and the exact 0{al) distribu¬ 
tions in Pt and Y for the qq reaction in the DIS scheme for the values of po chosen. However, 
for the gg channel in the MS scheme the resummation produces a large difference. The 
difference between the perturbation improved and the exact 0{al) distributions depends on 
the mass factorization scheme (DIS or MS), the factorization scale p, as well as the specific 
reaction under consideration {qq or gg). For a mass m = 175 GeV/c^ the gg channel is not as 
important numerically as the qq channel. However, since the corrections for the gg channel 
are quite large, resummation produces a non-negligible difference between the perturbation 
improved and the exact 0{al) distributions when adding the two channels. However, the 
shapes of the distributions are essentially unchanged. 

ACKNOWLEDGEMENTS 

We thank E. Laenen and W. L. van Neerven for helpful discussions. 

The work in this paper was supported in part under the contract NSF 93-09888. 


16 



REFERENCES 


[1] (CDF Collaboration) F. Abe et al., Phys. Rev. D 50, 2966 (1994). 

[2] (DO Collaboration) S. Abachi et al., FERMILAB-PUB-94-354-E, 1994. 

[3] L. H. Orr and W. J. Stirling, DTP-94-60, 1994. 

[4] L. H. Orr, T. Stelzer, and W. J. Stirling, DTP-94-112, 1994. 

[5] F. A. Berends, J. B. Tausk, and W. T. Giele, Phys. Rev. D 47, 2746 (1993). 

[6] R. H. Dalitz and G. R. Goldstein, Phys. Rev. D 45, 1531 (1992). 

[7] P. Nason, S. Dawson, and R. K. Ellis, Nucl. Phys. B303, 607 (1988). 

[8] W. Beenakker, H. Kuijf, W. L. van Neerven, and J. Smith, Phys. Rev. D 40, 54 (1989). 

[9] E. Laenen, J. Smith, and W. L. van Neerven, Phys. Lett. B 321, 254 (1994). 

[10] R. K. Ellis, Talk presented at the Glasgow Conference, July 1994. 

[11] A. D. Martin, R. G. Roberts, and W. J. Stirling, Phys. Lett. B 306, 145 (1993). 

[12] P. Nason, S. Dawson, and R. K. Ellis, Nucl. Phys. B327, 49 (1989); B335, 260(E) 
(1990). 

[13] W. Beenakker, W. L. van Neerven, R. Meng, G. A. Schuler, and J. Smith, Nucl. Phys. 
B351, 507 (1991). 

[14] E. Laenen, J. Smith, and W. L. van Neerven, Nucl. Phys. B369, 543 (1992). 

[15] K. Lane, BUHEP-95-2, 1995. 

[16] R. Meng, G. A. Schuler, J. Smith, and W. L. van Neerven, Nucl. Phys. B339, 325 
(1990). 

[17] G. Sterman, Nucl. Phys. B281, 310 (1987); D. Appel, P. Mackenzie, and G. Sterman, 
ibid. B309, 259 (1988); S. Gatani and L. Trentadue, ibid. B327, 323 (1989); B353, 183 


17 



(1991); H. Contopanagos and G. Sterman, ibid. B400, 211 (1993). 

[18] A. H. Mueller, Nucl. Phys. B250, 327 (1985); M. Beneke, Phys. Lett. B 307, 154 (1993); 
G. Grunberg, GPTH-A328-1094, 1994; Yu. A. Simonov, Pisma Zh. Eksp. Teor. Fiz. 57, 
513 (1993) [JETP Lett. 57, 525 (1993)]. 

[19] H. Gontopanagos and G. Sterman, Nucl. Phys. B419, 77 (1994). 

[20] L. Alvero and H. Gontopanagos, ITP-SB-94-33, 1994; H. Gontopanagos, L. Alvero, and 
G. Sterman, ITP-SB-94-40, 1994. 

[21] R. Bamberg, W. L. van Neerven, and T. Matsuura, Nucl. Phys. B359, 343 (1991); W. 
L. van Neerven and E. B. Zijlstra, Nucl. Phys. B382, 11 (1992); P. J. Rijken and W. 
L. van Neerven, INLO-PUB-14-94, 1994. 

[22] H. Plot how-Besch, ’PDFLIB: Nucleon, Pion and Photon Parton Density Functions and 
as Calculations’, Users’s Manual - Version 4.16, W5051 PDFLIB, 1994.01.11, GERN- 
PPE. 


18 



FIGURES 


(k) - 

FIG. 1. The top quark pt distributions da\j /dpt for the qq channel in the DIS scheme for a 
top quark mass m = 175 GeV/c^. Plotted are da^lP/dpt (upper solid line), da^^/dpt |exact (lower 
solid line), da^^ /dpt |app (upper dotted line), da^^ /dpt |app (lower dotted line), and da^ff/dpt 
{po = 0.05 m upper dashed line and po = 0.1 m lower dashed line). 

FIG. 2. The top quark pt distributions dan/dpt for the qq channel in the DIS scheme for a top 
quark mass m = 175 GeV/c^. Plotted are da^jp /dpt + da^^/dpt lexact (solid line) and da^^'^/dpt 
(po = 0.05 m upper dashed line and po = 0.1 m lower dashed line). 

FIG. 3. The top quark pt distributions dajj /dpt for the gg channel in the MS scheme for a top 
quark mass m = 175 GeV/c^. Plotted are da^^ /dpt (upper solid line at large pt), da^^/dpt lexact 
(lower solid line at large pt), da^^/dpt |app (lower dotted line), da^^ /dpt |app (upper dotted line), 
and da'^ff /dpt {po = 0.2 m upper dashed line and po = 0.25 m lower dashed line). 

FIG. 4. The top quark pt distributions dan/dpt for the gg channel in the MS scheme for a top 
quark mass m = 175 GeV/c^. Plotted are dafj^/dpt + da^^ /dpt lexact (solid line) and da^'^/dpt 
(po = 0.2 m upper dashed line and po = 0.25 m lower dashed line). 

FIG. 5. The top quark pt distributions da\j /dpt for the sum of the qq and gg channels for a top 
quark mass m = 175 GeV/c^. Plotted are da^^ /dpt (upper solid line), da^^ /dpt lexact (lower solid 
line), da^^/dpt |app (upper dotted line), da^^/dpt |app (lower dotted line), and da'^ff/dpt (upper 
and lower dashed lines). 

FIG. 6. The top quark pt distributions dan/dpt for the sum of the qq and gg channels for a top 
quark mass m = 175 GeV/c^. Plotted are da^jp/dpt + da^^/dpt lexact (solid line) and da^^^/dpt 
(upper and lower dashed lines). 
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(k) - 

FIG. 7. The top quark Y distributions da\j /dY for the qq channel in the DIS scheme for a 
top quark mass m = 175 GeV/c^. Plotted are da^^/dY (upper solid line), da^jp/dY [exact (lower 
solid line), da^^/dY |app (upper dotted line), da^^ /dY |app (lower dotted line), and da'^ff/dY 
{fiQ = 0.05 m upper dashed line and /Uq = 0.1 m lower dashed line). 

FIG. 8. The top quark Y distributions da^/dY for the qq channel in the DIS scheme for a top 
quark mass m = 175 GeV/c^. Plotted are da^j^/dY + da^jp/dY [exact (solid line) and da^^/dY 
{fiQ = 0.05 m upper dashed line and /Uq = 0.1 m lower dashed line). 

FIG. 9. The top quark Y distributions dajj /dY for the gg channel in the MS scheme for a top 
quark mass m = 175 GeV/c^. Plotted are da^pdY (upper solid line at T = 0), dappdY [exact 
(lower solid line at T = 0), dappdY [app (lower dotted line), dappdY [app (upper dotted line), 
and da^ff /dY (^o = 0.2 m upper dashed line and go = 0.25 m lower dashed line). 

FIG. 10. The top quark Y distributions dau/dY for the gg channel in the MS scheme for a top 
quark mass m = 175 GeV/c^. Plotted are da^pdY + dappdY [exact (solid line) and da^^/dY 
{go = 0.2 m upper dashed line and go = 0.25 m lower dashed line). 

FIG. 11. The top quark Y distributions dajj /dY for the sum of the qq and gg channels for a 
top quark mass m = 175 GeV/c^. Plotted are dapp/dY (upper solid line), da^pdY [exact (lower 
solid line), dappdY [app (upper dotted line), dappdY [app (lower dotted line), and da'^ff/dY 
(upper and lower dashed lines). 

FIG. 12. The top quark Y distributions dan/dY for the sum of the qq and gg channels for a 
top quark mass m = 175 GeV/c^. Plotted are dappdY + dappdY [exact (solid line) and da^^^/dY 
(upper and lower dashed lines). 
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